Crystalline symmetry and the melt-growth kinetics of solid-liquid interface 
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One of the important factors governing the growth morphology of materials is the interface kinetic 
coefficient /i, which is the proportionality constant between the velocity of solid-liquid interface and 
undercooling. We employ Ginzburg-Landau (GL) free energy functional to derive an analytical 
expression of kinetic coefficients. The anisotropy of kinetic coefficients naturally arise from the 
broken symmetry at the solid-liquid interface for various crystalline orientations. The analytical 
expression of kinetic coefficients is compared to Mikheev-Chernov theory [J. Cryst. Growth 112, 
591 (1991)] derived from hydrodynamic equations. In addition, we use equilibrium density wave 
profiles to evaluate kinetic coefficients and compare them with that from MD simulations. Our 
results are in good agreement with Mikheev-Chernov theory and MD simulations and shed lights 
on a possible origin of anisotropy of interfacial kinetics. 
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I. INTRODUCTION 

The growth of crystal phases from the melt occurs in 
many material processing applications such as zone re- 
finement and epitaxial growth. In the early 20th cen- 
tury, a classical theory by Wilson for the growth of crys- 
tals rate from viscuous melt is proposed based on the 
picture of abrupt interface separating semi-infinite melt 
and crystal pQ. It assumes that the particles in the liq- 
uid phase diffuse over a distance approximately of the 
order of the particle size before impinging on and later 
being absorbed by the lattice. Comparing to experiments 
on crystal growth from viscous melts, this theory shows 
good agreements in terms of the crystal growth ratejSJ- 
7 . However, recent studies on crystal- melt interface of 
organic crystals p], simple metals [3], and Lennard- Jones 
crystals 4J have provided more insights into the growth 
kinetics and morphology of materials. In contrast to the 
growth from viscous melts, which consists of large com- 
plex molecules, the above-mentioned crystal-melt inter- 
face grows at a much higher rate than that predicted by 
the classical theory [21 [3]. Mikheev and Chernov |131 RB] 
first offered an analytic treatment of this interface growth 
kinetic by taking proper account of the crystal ordering 
at the crystal-melt surface. They propose that as the 
crystal grows the ordering of atoms occurs simultane- 
ously in many surface layers instead of solely on a sharp 
crystal-melt interface. With this notion, they derive an 
analytical expression for kinetic coefficient, which is the 
proportionality constant between the crystal-melt inter- 
face velocity and the interface undercooling, in terms of 
bulk physical quantities such as relaxation time, struc- 
ture factor, latent heat of fushion and melting tempera- 
ture. 

In the past few decades, several microscopic theories 
have been proposed to show that anisotropy of the excess 
interface energy is a crucial parameter that determines 
morphology of dendrites during solidification. Ginzburg- 
Landau (GL)theory with order parameter representing 
amplitude of density waves has been used to calculate 



the magnitude of interfacial free energy 7 in different sys- 
tems and also for different crystal structure like fee and 
diamond cubic [TTl [30] • GL theory is shown to suc- 
cessfully capture capillary anisotropy of crystal-melt in- 
terfaces by incorporating underlying crystal symmetries. 
In this paper, we derive interfacial kinetics using GL 
theory from a different yet more compact perspective. 
Firstly we outline the hydrodynamic approach adopted 
by Mikheev and Chernov [131 [H] in deriving the kinetic 
coefficient. Then we extend Ginzburg-Landau theory to 
describe a moving interface of solid-liquid systems un- 
der a given undercooling. Based on our previous work 
on crystalline anisotropy for bcc-liquid interface, we ex- 
plicitly computed the kinetic coefficients using the equi- 
librium density wave profiles, which is contrasted with 
that obtained through the approximation scheme pro- 
posed previously [121 Q3]. Finally, Comparisons to molec- 
ular dynamics(MD) simulation and the physical origin of 
kinetic anisotropy are discussed. 

Outlines of the Hydrodynamic approach: To under- 
stand the kinetics of the non-equilibrium solid-liquid in- 
terfaces of a given undercooling, it is instructive to con- 
sider the number density current j(r,t) (i.e. number of 
particles flowing across an unit area in a given direction 
per unit time) as well as the quadratic isothermal dissi- 
pative function [XT] IT3] constructed accordingly: 



R= d 3 rd 3 r'T a ^\{n})j a (r)j^), (1) 



where T a p is the (generally nonlocal and density depen- 
dent) dissipated kernel. The dissipated function R rep- 
resents the amount of (free) energy dissipated per unit 
time. Assuming that the heat transfer at the relevant 
length scale is much faster than the structure rearrange- 
ment under consideration !22, 23 j: the generalized hydro- 
dynamic equations are|llj: 
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Note that F[{n(r)}] refers to the free energy for the (in- 
homogeneous) equilibrium fluid minus its kinetic energy 
contribution^^, which is different from the free energy 
in the GL theory used later. 

For small undercoolings, the density profile for the 
solid-liquid interface can be described by the contant ve- 
locity wave of the form |13l |2"3] : 
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where the sum is taken over principal reciprocal lattice 
vectors of the crystal. Note that we align the inter- 
face normal in z, and denote the steady-state velocity 
as V. By solving Eqs.QQQ one immediately gets 
the free energy balance relation between dissipation and 
crystallization 13J : 
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with the dissipated kernel given by 
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. Here no represents the density of liquid, T m the melt- 
ing temperature, and AT = T — T m the undercooling 
with reference to the melt. Note that the dissipated ker- 
nel is related to the relaxation time T(Ki) measured as 
the inverse half-width of the dynamical structure factor- 
s'^, w) at given k = Kj\2A . Inserting the longitudinal 
part of the traveling density wave Eq. Q into the energy 
balance equation Eq. ^ one arrives at the expression for 
interface velocity: 
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As will be shown latter, Eq.([7]) can be derived from the 
GL theory by including a linear temperature dependence. 

Extension of Ginzburg- Landau Theory: Different from 
the Ginzburg-Landau theory for the solid-liquid interface 
at equilibrium developed by Shih et al.\18\ and Wu et 
aL|17j. we extend the discussion to a solid- liquid system 
that is out of equilibrium where the system is undercooled 
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FIG. 1: (Color online) Density wave profile Eq. Q for inter- 
face moving to the right with velocity V . 



and the solid-liquid interface moves at a constant veloc- 
ity. Near the melting point T m , the free energy difference 
between solid and liquid phase is related to latent heat 
of fusion L, 



F S {T) - F L (T) = {S s - S L )(T - T m ) = L 
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T 
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where we have used the thermodynamic relation dF = 
—SdT. Since we are to deal with the dynamics close to 
the melting temperature, we should incorporate a term 
linear in u (i.e. the amplitude of density waves corre- 
sponding to a specific principal lattice vector of the re- 
ciprocal lattice of the solid) with the equilibrium free 
energy so as to recover Eq,(|8|). Thus, one can write 



AF' = AF- 
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where AF is the excess free energy for the solid-liquid 
system in equilibrium, which, under the assumption that 
the density wave amplitude varies slowly over a scale 
I /F m ax where K max represents the wave vector at which 
the peak of structure factor is located, takes the form[T7]: 
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Note that C(K) refers to the Fourier transform of the 
direct correlation function C{\r — f|) and C"(K) = 
d 2 C(K) /dK 2 , n Q is the atomic density in the liquid 
phase, and u s is the equilibrium amplitude of the prin- 
cipal reciprocal lattice vectors in the solid which corre- 
sponds to one minimum of AF(u,T m ). One can check 
that this is point where AF' (u s ,T m ) = 0. Note that 
for convenience in subsequent comparisons we include a 
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normalization factor 1/6 to account for the slope of the 
linear tilt potential of the bcc-liquid interface [17], 

As usual, one can obtain the kinetic coefficient v by 
looking for the steady-state propagating solutions of the 
following equation of motion(assuming the interface nor- 
mal n = z) 

du t dut SAF' 

where r = r(G) / S(G) is the relaxation time mea- 
sured as the inverse half-width of the dynamic structure 
factor [21]. An analytic expression for V and hence /x, the 
kinetic coefficient, can be obtained by using a solvability 
condition as outlined in all phase-field review papers J25J- 



|2"7] . which relates V and T — T m . For anisotropic cases, 
the interface profiles for different sets of density wave am- 
plitudes are coupled but /j, can still be computed analyti- 
cally from a similar solvability condition whose procedure 
is shown below. 

Our previous work on crystalline anisotropy for bcc- 
liquid interfaces [17] shows that the amplitudes of density 
waves can be categorized into different groups according 
to the value of square gradient term appeared in the GL 
free energy functional(summarized in TABLE. I). For this 
reason, we arbitrarily choose the {110} crystal face along 
which the amplitude of propagating density waves are 
denoted as u,v, and w to examplify our calculation. 



Different orientation of crystal faces 


100 


110 


111 


(Ki-n)\ 
Number of K[s 
Square gradient term c, = (Ki ■ n) 2 /A 


1/2 
4 8 
1/8 


1/4 1 
8 2 2 
1/16 1/4 


2/3 
6 6 
1/6 



TABLE I: Classifications and the values of square gradie nt a of different orientations of bec crystal interfaces. 



Using the equation of motions Eq.( 11 ) we have 
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where a is the linear tilt appended to the the free energy 
of the equilibrium case 
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To obtain a steady-state solution of propagating wave 
with velocity V along the interface normal z, we cast 
Eq.(12| into the moving frame and expand the density 



wave amplitude around the melting point. After using 
the stationarity of the equilibrium density and simplify- 
ing, one arrives at the following matrix equation: 
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Here we denote uq, vq and wq as the equilibrium density 
wave amplitudes specific to the {110} interface. Note 
that AF — J drf(u,v,w), f u = Sf/Su. Employing again 
the stationariy of equilibrium density profile, we have 
£ty = o with V = (^,^,^) T - Combining this 



with Eq. ( 14 ) one arrives at 
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Setting the boundary conditions uo(— oo) = vq(— oo) = 
Wq(— oo) = and Mo(oo) = v (co) — wq(co) — u s , the 
density wave velocity V can be further simplified into 
(here the subscript of V indicates the crystal face normal 
{110} specific to this case) 
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FIG. 2: (Color online) Velocity versus dimensionless under- 
cooling AT = (T — T m )L/fcsT^ l for different crystal faces 



We can also apply the same analysis to different crys- 
tal faces with references to different sets of density wave 
amplitudes and square gradients listed in TABLE. I. In 
any case, we can re-write the kinetic coefficient (/z) in 
a more compact form by noting that V — /iAT and 
AT = (T — T m ): 



As one can see, this is the same as that derived by 
Mikheev et'al [T3| (see Eq.Q). We plotted the depen- 
dence of velocity V with reference to different crystal 
faces on the dimensionless undercooling AT = (T — 
TjnjL/kBT^ explicitly in FIG(2]and contrasted the nu- 
merical value of \x obtained here with that in MC theory 
in TABLE II. Although both theories predict similar ki- 
netic anisotropy ordering[28j. doo > Mno > Mill) we 
would like to emphasize that MC theory (i.e. the sec- 
ond row in TABLE II) computed the kinetic coefficient 
with Ui estimated from a linear theory, which is differ- 
ent from the exact profiles we used here(i.e. the first 
column in the same table). The good agreement of our 
calculation with both MD and MC theory reveals that 
by taking into consideration the anisotropy of crystalline 
orientations the broken symmetry of interfacial kinetics 
is automatic regardless of (i)the details of atomistic dy- 
namics within or outside of the interface as well as (ii) 
the hydrodynamcal properties of the melt. 
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GL Theory 


Mikheev-Chernov 


MD 


Mioo/mho 


1.06 


1.14 


1.37 ±0.07 


Mioo/ Cm 


1.27 


1.23 


1.44 ±0.07 



TABLE II: The anisotropy of the kinetic coefficient for bec lattices calculated by GL theory derived in this work, 
Mikheev-Chernov model [13], and MD simulation with MH( SA) potential [Wj IT?) 
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Conclusions: We have computed the kinetic coeffi- 
cient /x, which is the proportionality constant between 
interface velocity and undercooling, by incorporating a 
temperature dependence term into the Ginzburg-Landau 
(GL) free energy functional. Though different from the 
Mikheev-Chernov (MC)[T3J theory based on the hydro- 
dynamic equations of conservation of momentum and 
particle number, our intuitive and more compact exten- 
sion of Ginzburg-Landau calculation predicts the same 



ordering for kinetic coefficients in different crystal faces. 
In addition, we have contrasted both theories with MD 
simulations and find good accuracy. This reveals that 
not only is our method quantitatively precise in dealing 
with non-equilibrium interfacial kinetics problems but 
also provides a explanation to the anisotropy of inter- 
facial kinetics - the broken symmetry at the interface for 
various crystalline orientations. 
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